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2+1 flavour DWF simulations 



Peter Boyle 



The RBC and UKQCD collaborations have jointly performed 2+1 flavour simulations of QCD, 
representing the up, down, and strange quarks with the standard domain wall fermion action. This 
review will firstly cover the theoretical foundations of our simulations considering issues such as 
locality, chirality and topology. We will secondly present the ensemble parameters and algorithms 
used in our simulations, and discuss algorithmic performance and trade-offs. We will finally sum- 
marise important results presented elsewhere in this conference [|l|, |2[ Isj, ||, ^ and recent publi- 
cations [11, 12. I^, 0' iJ' 0' iZ' 1^' l^l- These results include pseudoscalar masses and decay 
constants and low energy constants of the chiral effective lagrangian. We also review results for the 
neutral kaon mixing amplitude Bk, the KB form factor, pseudoscalar meson structure, and vector 
meson decay constants. In the baryon sector we review results for the spectrum, and nucleon form 
factors and structure functions. Highlights of our programme include preliminary quark masses, 
and determinations of V„.s from both fx/ fn and from K13, and an updated result for Bk- We find 
significant finite volume effects in the nucleon axial charge gA for our mj^ = 330 MeV ensemble on 
a (2.7fm)^ lattice, and highlight the importance of large physical volumes for non-trivial nucleon 
physics. 



1. Lattice action, algorithms, cost 

We use the Iwasaki gauge action and domain wall fermion action 

dI.{m,)=d^{-Ms) 



DtArnf) = 2[(l-75)5.+iy + (l + 75)5.-iy-25,^ 



2 



(1 - 75)5s-,L,-l^,y + (1 + 75)5.,o5l,-1/ 



(1.1) 



Here Dyv is the Wilson Dirac operator, and the boundary conditions are understood to be Dirichlet 
in the fifth dimension, periodic in spatial directions and anti -periodic in time. Surface states of 
either chirality are bound to the 4-dimensional s = and s = Ls — I hyperplanes and are identifed 
with physical, four-dimensional modes 

qix) = Pl»P(^,0) +Pr'¥{x,L, - 1). 

In dynamical simulations the bulk infinity of the five dimensional partition function is removed 
using Pauli-Villars fields. 

Our simulations, table [l], have principally been performed using a single value for the lattice 
spacing, a^' = 1.73 GeV, and using both 16^ and 24^ lattice volumes corresponding to (2.0fm)^ 
and (2.7fm)^. A status report on simulations in progress at a second, finer lattice spacing with a 
32^ lattice is also given. These latter simulations are carried out as part of a collaboration between 
RBC, UKQCD, and LHPC. 

The degree of flexibility and choice of the implementation of Hybrid Monte-Carlo has ad- 
vanced greatly in recent years, with several new algorithmic variants proposed. These advances 
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L^xTx L, 


{ami,ms) 




a-^ (GeV) 


L(fm) 


nin (MeV) 


"^res 


tMD 




(0.01,0.04) 








400 




4000 


16^ X 32 X 16 


(0.02,0.04) 


2.13 


1.62(4) 


1.94 


530 


3.08 X 10-^ 


4000 




(0.03,0.04) 




(P) 




630 




7500 




(0.005,0.04) 








330 




4500 


24^ X 64 X 16 


(0.01,0.04) 


2.13 


1.73(3) 


2.73 


420 


3.15 X 10-3 


4700 




(0.02,0.04) 




(Q-) 




560 




2800 




(0.03,0.04) 








670 




2800 


32^ X 64 X 16 


(0.004,0.03) 


2.25 


~2.15 


~2.93 


~260 


~ 6 X 10-4 


1100+ 




(0.006,0.03) 








~310 




1300+ 



Table 1: Ensemble parameters for the UKQCD/RBC data set. The 32^ ensemble production is in collabo- 
ration with LHPC since July 2007. The negative Wilson mass in the domain wall formalism was 1.8 for all 
ensembles. 



include the (affordable) extension of exact algorithms to odd numbers of flavours [|20||, several 
schemes for splitting the fermionic force into UV and IR portions that can be updated on different 
timescales, and improved numerical integrators. Experimenting with the available options^ we 
have settled on RHMC with a hybrid combination multi-mass preconditioning at light mass scales 



[22] and multiple pseudofermion fields at [21] heavier mass scales with Omelyan integrators [L 



The simulated fermion determinant is included as 



D'^D{m,)] B \ D(l) ] B \ D(l) ] B \ D(l 

where each determinant factor is estimated via a separate pseudo-fermion field (thus four in all). 
The degenerate u,d flavours are mass preconditioned by the strange mass, the remaining three 
factors of the strange mass make use of the RHMC n-roots force reduction trick, and the factors of 
D(l) are the Pauli-Villars fields. The "A" and "B" determinants are updated on different timescales 
using nested Omelyan integrators with Omelyan parameter A = 0.22. 

The coarsest timescale used is 5^ = g for the most expensive up/down fields at our lightest 
mass. A trajectory length T = 1 is used and thus contains only six timesteps. The nature of the finer 
timesteps are somewhat complicated by our use of the Omelyan integrator. An Omelyan integra- 
tion QPQPQ timestep involves two force calculations that are not equally spaced in Monte-Carlo 
time. Reversibility leaves little flexibility for possible approaches to integrator nesting and in our 
nomenclature a 1 : 1 nesting implements a complete QPQPQ Omelyan timestep of the second force 
contribution for each sub-timestep of the first. Our strange mass determinants are 1:1 nested inside 
the up/down determinants implying that each of the three "B" force contributions are calculated 
twice for every "A" force contribution. The gauge force is nested inside the "B" contribution in a 
similar way but with an Omelyan nesting ratio of 1:6. Convergence residuals used in molecular 
dynamics phases vary between 10^** and 10^^ according to the typical force contribution, while 
lO^'" is used uniformly for all Metropolis steps. Guesses are history independent and reversibility 
has been demonstrated to very high precision. 



we gratefully acknowledge the immense contribution of Mike Clark to our programme 
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The two most important measures of cost are the technology independent algorithmic cost and 
the wall clock time to run an ensemble on the machines available to the collaboration. The scaUng 
with light quark mass is weak for our mass-preconditioned algorithm. To allow concrete com- 
parisons with other calculations, we quote that the 24^ ensemble with m/ = 0.02 requires (9(10^) 
applications of for a T = 1.0 trajectory. Around 40 T = Itrajectories per day are produced on 
a 4096 node QCDOC partition and sustained performance is around l.lTFlop/s on this machine 
size. For our 32^ simulations on the same machine size around 10 units of MD time are produced 
per day, with a trajectory length of t = 2 MD time units. 



2. Theoretical foundations 



The DWF five dimensional system can be represented as a Fock space trace with a transfer 
matrix T = e~^'^ where tanh ^ = 2+i)w ~ ^^'^ dimensional effective action of DWF is a 
functional of the gauge fields and is not manifestly local. This approximates an overlap operator 
making use of Hj as the argument to the sign function. 

[det£)dwf(l)]~^det£)dwf(m) = det^ [l +m + 75(l -m)tanh(L^tanh"^*:5)] 

^ det^[l+m + (l-m)75Sgnii:s] (2.1) 

Ignoring anomalous chiral symmetry breaking for now, consider the flavour non-singlet axial 
current in this formulation. The five dimensional theory has a conserved five dimensional vector 
current. The DWF (five dimensional) axial transformation associates positive and negative chiral 
charges with the positive and negative halves of the fifth dimension. One can construct a four 
dimensional axial current that is extensive in the fifth dimension and for which the chiral symmetry 
breaking effect of finite Lg consists only of a mid-point term in the fifth dimension. 

where 

<W= l'sign(.-^)yj(x,.), 

/"'(x) = ^(x,0)PRt'''¥{x,L, - 1) - »F(x,L, - l)/i?"»P(x,0) = q{x)Y5t''qix) 

y|,(x) =^(x,L,/2)Pr?Xx,L,/2- 1) -^(x,L,/2- l)Pi?Xx,L,/2) 

In low energy Greens functions, the midpoint density J^^ is equivalent to the dimension-three op- 
erator 2mres^y^^ where mres is an additive mass renormahzation, measured as 

{Jl^{x)ny)) 



2(P«(x)P°(y))' 



All unphysical chiral symmetry breaking effects in DWF, including mres, involve propagation from 
a source field of one or other chirality on the corresponding domain wall across the fifth dimension. 
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These naturally involve the transfer matrix T raised to an appropriate (large) power. This suppres- 
sion mechanism is key to the quality of DWF lattice chiral symmetry, and thus understanding the 
details of the nature of the spectrum of Hj are crucial. The asymptotic propagation in the fifth 
dimension is, for large L^, dominated by the modes of Hj at the lowest eigenvalues at which the 
eigenmode density p(A) is non-zero. However, for modest Lg these asymptotic contributions can 
be very much suppressed and details of Ht and the size distribution of the eigenmodes must be 
considered. 

The translational invariance of the DWF approach in the fictitious fifth dimension admits 
power counting in T^' as a powerful tool, and wires serves as qualitative guide to the cost of one 
factor of T^'. It is worth emphasizing that this counting carries real power; for example wrong 
chirality mixings for Bk require two crossings of the fifth dimension and are proportional to m^^^ 
[^7|]. Calculations of |- are feasible [28, 8|], while residual chiral symmetry breaking is sufficiently 



enhanced in direct measurement of the chiral condensate that it will prove problematic with DWF 
without further reduction in rrires (or matching large Ls valence simulations to our sea pion masses) 



[28] 



2.1 Localisation structure of Aoki phase 

The zero modes of Hj and Hw coincide and the known quenched Aoki phase behaviour of 
H]/v is directly relevant. It has recently been understood that the structure of the quenched Aoki 
phase is rich: there is a non-zero density of near zero modes of H^, and a Banks Casher pionic 
condensate, throughout all of the accessible phase diagram; however, a localisation transition is 
thought to occur and, towards weak coupling, the phase is thought to display a non-zero mobility 
edge. Establishing that we at least have a non-zero mobility edge (or better yet a gap in the spec- 
trum) is key to establishing locality of dynamical overlap or DWF simulations. The conjectured 
model for the structure of the spectrum has p(0) 7^ and p(Ac) ^ p(0) - a volume factor enhances 
the contributions from modes above the mobility edge A^. A consequence of this model for the 
structure of the spectrum of Hj is that 25] 



Ls 

Here these two contributions come from a low density of (volume factor suppressed) localised 
near-zero modes and a larger density of extended modes near the mobility edge. The overall factor 
of 7- represents an infrared cut off on the shell of modes that contribute significantly imposed by 
the size of the fifth dimension. 

2.2 Locality 

For sufficiently smooth gauge fields there is a gap in the spectrum of and this implies 
locality of the corresponding overlap Dirac operator taking H\y as the kernel of the sign function 
[p4|]. This proof may be generalised to cover the Shamir Kernel that corresponds to the — > 00 
limit of DWF, and also the condition can be relaxed to require only a gap in the spectrum of 



delocalised eigenmodes [g^, |13[ ]. 

As there is no gap in the spectrum of this Hj (in the absence of a ghost Wilson determinant) it 
is necessary to demonstrate that its spectrum displays a non-zero mobility edge to establish locality 
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L 

Figure 1: We display the dependence of niies on Ly for valence quarks on Domain Wall 2+ If ensembles with 
the Iwasaki gauge action and — 8. The fit correspond to a model based on the conjectured mobility edge 
structure in the Aoki phase, and a good description of our data is seen. This evidence of a non-zero mobility 
edge indicates locality for the DWF effective action. 

of the effective four dimensional theory. The related matrix = YsD^ which has identical zero 
mode structure, and it also suffices to study Hw in its place. 

This can be done in two ways. Firstly we have done so indirectly by demonstrating the consis- 
tency of the behaviour of wires (Av) with the above model, figure |l} Secondly microscopic inspection 
of the eigenmodes of Hw can check the locality of individual low lying eigenmodes on a mode by 
mode basis, figure ^. Given an eigenmode, we take y as the location of the maximum 

of Y^Yi^) ^rid firid the lowest exponential localisation length Leff that for all coordinates x with 
\^ — y\ ^5 satisfies bound 

2|.v-y| 

We note that the Leff bounds the eigenmode large distance from its peak in all directions, and 
emphasize that this strict bound approach is robust against eigenmodes with extended lower di- 
mensional sub-spaces and other pathological cases. We therefore have demonstrated a non-zero 
mobility edge Xc > 0.2, and thus locality, for our jS = 2. 13 simulations and that we are therefore in 
a correct part of the Aoki phase diagram for taking a continuum limit. 

2.3 Chiral symmetry breaking and non-perturbative renormalisation 

We use the Rome-Southampton Rl-mom approach to determine the renormalisation of our 
lattice operators non-perturbatively. The valence DWF action suppresses 0(a) effects, both on and 
off shell, and thus is particularly well suited to the off-shell renormalisation approach. The good 
chiral symmetry of DWF is reflected in its renormalisation structure, and one might expect this to 
be well demonstrated and tested by NPR. 

In practice, however, there are substantial physical spontaneous chiral symmetry breaking 
effects at low for standard Rl-mom kinematics. These obscure a demonstration of the good 
chiral properties of DWF in coarse lattice spacing simulations. An example of this for the symmetry 
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Figure 2: We look at the scatter of exponential bounds measured for low modes of the Hermitian Wilson 
Dirac operator forMs = — 1.8 for j3 = 2.13. A mobility edge is clearly seen, giving direct evidence of the 
locality of DWF for our simulated parameters. 




Figure 3: The left panel shows a class of vertex correction graphs admitting non-perturbative physics accom- 
panied by only a^/ suppression. This is a disadvantageous feature of the standard Rl-mom kinematics. 
The right panel displays the (physical spontaneous) chiral symmetry breaking effects that split the axial and 
vector amputated vertices — Ay for these kinematics. The splitting is correspondingly poorly suppressed 
in l/p^. 
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MS 2 GeV 


Z 


Stat 


sys 


Za 


0.7161 


0.0001 




7 


1.656 


0.048 


0.104 


Zt 


0.7951 


0.0034 


0.0117 




0.9276 


0.0052 


0.0222 



Table 2: Renormalisation constants for our /3 =2.13 ensembles. Results are quoted in the chiral limit, and 
obtained on the 16^ ensembles. Systematic errors quoted contain estimates of extrapolation errors and those 
from (continuum) perturbative conversion to MS. The exception is which is an improved measurement of 
the ratio of the conserved axial current ^ (x) to the more commonly measured boundary field axial current 



breaking splitting of — Ky in shown in figure |[ These effects also introduce an ambiguity in the 
determination of the renormalisation constants of around 2%. This physical effect is three orders 
of magnitude larger than any contamination expected from our residual chiral symmetry breaking 
"^res = (3 X 10^^)^ , and no improvement will be gained by improving the accuracy of the chiral 
symmetry by either increasing or using a more exact overlap approach. 

The problem has been enhanced by the particular choice of kinematics used in traditional RI- 
mom NPR. For a standard bilinear vertex function the leg momenta are equal and a soft subgraph 
is only suppressed by the of a single hard gluon as shown in figure ^. 

A better alternative is to gain further suppression of soft contributions using non-exceptional 
momentum kinematics = p'^ = {p — p'Y, figure^ 

Unfortunately the large body of higher perturbative calculations for various operators in RI- 
mom is non-trivial to reproduce for these alternative kinematics. However, use can be made of these 
results since even without conversion functions the good chiral chiral properties of DWF can be 
demonstrated and omission of chirality mixing for B/f justified without requiring the perturbative 
conversion. There is no practical benefit, for these specific quantities, from further redution in 
mies- It is better to focus our available effort on more pressing problems. The great improvement 
in the scaling window we have demonstrated for Rl-mom vertex functions with non-exceptional 
kinematics should serve as encouragement to the revisit high order calculations of the conversion 
to MS for non-exceptional kinematics. It is also quite possible that the convergence properties 
of the perturbative expansion will be improved in some cases. For example the four-loop Wilson 
coefficient relevant to displays much better convergence in MS than for conventional Rl-mom. 

Table § lists the axial cun^ent, field, mass, tensor and four quark operator renormalisation 



constants obtained in reference |38]. 



2.4 Topological tunneling 

One of the principal attractions of dynamical fermion simulations with good chiral symmetry 
is the existence of a correct axial anomaly. However, this perfection of the action could easily be 
compromised by mundane algorithmic issues resulting in failure to sample the topological structure 
of the vacuum adequately. Problems arise with exact chiral symmetry since the sign function 
has a discontinuity which molecular dynamics updates will skip over for any non-zero timestep. 
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Figure 4: By using non-exceptional momenta, spontaneous chiral symmetry breaking effects are pushed 
to much lower momenta than with standard Rl-mom kinematics. The right panel displays a cross-chirality 
mixing matrix elements relevant to Bk- The excellent chiral properties of DWF, when extrapolated to the 
chiral limit are now apparent even at intermediate momentum scales. This approach is promising and can 
demonstrate the absence of unwanted lattice mixings. In order to reduce NPR systematics for renormaUsed 
quantities recalculation of Rl-mom anomalous dimensions is required to high order in continuum perturba- 
tion theory for these new kinematics. The left panel displays the difference between the amputated axial and 
vector vertex functions. 



Exact overlap 
DWFL^=16 



-0.5 



Figure 5: Tanh approximation to the sign function forL, = 16 



Approximations to the sign function, DWF included, can involve a smooth transition over some 
eigevalue range, figure ^. Problems of this nature have not been observed with DWF simulations, 
and are unlikely to appear as the integration is problem free provided dz <^ Healthy global 
topological charge histograms are obtained on our 24^ ensembles in figure ^ and the susceptibility 
is displayed in figure 

The molecular dynamics problems have resulted in two responses from the dynamical overlap 
community. The reflection-algorithm treats the lowest modes exactly using a particularly expensive 
approach, while others [^] have used an auxiliary pseudofermion determinant to freeze the global 
topological charge by suppressing the density of low modes of Hw- Here, it is worth noting that the 
Zolotarev approximation to the sign function used has a coarse lower bound | A | > 0. 1 and does not 
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differ greatly from the tanh displayed in figure ^ In the language of the overlap, the key difference 
from the RBC-UKQCD approaches is that the authors of [29] suppress the density of low modes of 
Hw in the region —0.1 < A < 0.1 using an auxiliary fermion determinant, and it is then feasible to 
project and treat exactly the few remaining low eigenmodes of H^. This results in a fixed topology 
simulation with improved chiral symmetry. Similar approaches could equally well be combined 
with the DWF tanh approximation [||, ||, |3|, |4|]. 

We estimate that DWF simulations are around five times cheaper than the five dimensional 
Zolotarev approach used in [29], and twenty times cheaper than the nested four dimensional ap- 
proach. 

As DWF shows, issues with integrability of the fermion contributions in molecular dynamics 
are likely either solvable or avoidable. We note that were the auxiliary determinant and projection 
dropped, then the remaining Zolotarev approximation would yield a very "domain-wall-like" sim- 
ulation and tunneling would likely take place. Some thought has recently been given to developing 
algorithms that continue to tunnel topology in the presence of an auxiliary determinant [30]. 

A more fundamental problem is the increasing potential barrier introduced between topologi- 
cal sectors by the gauge action with increasing j3. Eq. 2A and discussion indicates that m^f^^{Ls = 
32) can be a qualitative guide to the near-zero mode density p(0). Figure^ shows that this, and we 
conjecture the tunneling rate per unit lattice volume, vanishes exponentially in the gauge coupling. 

The trade off is clear: suppressing the low modes that mediate topology change and attempting 
to answer the difficult question about ergodicity of the simulation is one rational choice. RBC and 
UKQCD's choice has been to accept a level of residual chiral symmetry breaking that is merely 
a minor irritant, but avoid questionable ergodicity and the risk of getting anomaly physics wrong 
in a particularly expensive fashion. This is also a substantially cheaper simulation, allowing more 
important systematic issues to be addressed. 



3. Measurements 

Table |3| gives a summary table of the mesonic measurements made on our ensembles. Further 
measurements have been made of nucleon two and three point functions, and static-light two and 



three point functions [^, [T^]. Between 150 and 700 measurements have been made on each en- 
semble depending on the quantity, and valence pseudoscalar masses vary between 240 MeV and 
750 MeV. Dynamical pion masses run from 330MeV to 650 MeV. Many more valence masses 
than dynamical masses are used to increase the amount of information in the chiral regime and are 



exploited in fits to partially quenched chiral perturbation theory [|36|], as highlighted by the masses 
quoted in red. 



4. Chiral effective lagrangian 

In this section we review the results presented in more detail by Meifeng Lin and Enno Scholz 
at this conference [Q]. Two approaches to fitting our data for pion and kaon masses and decay 
constants to obtain the EEC's of the chiral effective Lagrangian have been presented at this con- 
ference. The finite range of validity of chiral perturbation theory leaves such determinations from 
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m =0.005 



m =0.01 



0.07 r 




Figure 6: Topological charge distribution on our four 24^ ensembles. These are m/ = 0.005 (top-left), 0.01 
(top-right), 0.02 (bottom-left), and 0.03 (bottom-right). The light fermion mass is clearly constraining the 
distribution, and we are Ukely sampling topology well enough to reproduce 0=0 QCD. 



lattice (or indeed real world) data as something of an art. The rather massive real world kaon is 
neither unambiguously light nor heavy compared with chiral scales. 

One approach is to fit the full SU(3)xSU(3) chiral effective theory to data including the kaon 
as an active chiral pseudoscalar. A strength of lattice formulations, such as DWF, with controlled 
flavour symmetry is that the chiral perturbation theory can be decoupled from lattice artefacts in 
our simulated flavour content. 

A good available alternative is to treat the kaon as a non-Goldstone boson, coupled to an 
effective SU(2) x SU(2) theory. The analysis is applicable whenever m„ <^ mx, whether or not the 
kaon is heavy or light compared to other scales and we do not rely on chiral perturbation theory 
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Iwasaki gauge action, a~0.13 fm 




Figure 7: Left panel displays topological susceptibility on our 24^ and 16^ configurations; the line is the 
leading order chiral behaviour with E taken from the Gell-Mann-Oakes-Renner relation. Right panel displays 
the dependence of mres(^i = 32) as a function of inverse gauge coupling. As discussed in the text this can be 
taken as loosely indicative of trends in the density of near zero modes and also of the topological tunneling 
rate per unit lattice volume. 



mx,my 


0.001 


0.005 


0.01 


0.02 


0.03 


0.04 


0.001 


A,B ; 240 


290 


340 


420 


490 




550 


0.005 


A,B 


A,B,C; 330 


370 


450 


520 




580 


0.01 


A,B 


A,B 


A,B,C;410 


480 


560 




600 


0.02 


A,B 


A,B 


A,B 


A,B,t;550 


600 




650 


0.03 


A,B 


A,B 


A,B 


A,B 


A,B,C; 650 




700 


0.04 


A,B 


A,B,C 


A,B,C 


A,B,C 


A,B,C 




A,B,C; 750 



Table 3: We display the available mesonic measurements for the four 24^ ensembles m/ G 
{0.005, 0.01,0. 02, 0.03};oti = 0.04. The table lists valence masses rux and niy. In the upper right trian- 
gle we give the approximate pseudoscalar meson mass in MeV composed of quarks with valence masses 
and my. Those masses quoted in red survive our cuts for making NLO partially quenched chiral fits. In the 
lower left triangle we denote meson mass and decay constant measurements by "A", neutral meson mixing 
matrix element measurements by "B", and semileptonic decay matrix elements and distribution amplitudes 
by "C". Black corresponds to valence measurements made only on the Ughtest two ensembles. Blue corre- 
sponds to valence measurements made on all four ensembles. Green corresponds to unitary measurements 
made only with the valence quark masses equal to sea quark masses. 

being convergent at kaon masses. The LEC's will be strange mass dependent and, since the kaon 
is somewhat lighter than a typical chiral scale, the convergence of the chiral expansion controlled 
by these LEC's may be correspondingly impacted. This will merely reflect the new dynamics that 
enters at the kaon mass scale. A reasonable estimate is that successive orders in the chiral expansion 
will only be suppressed by mi/ms. To estimate possible systematic NNLO contamination entering 
when we perform SU(2) fits with a mass cut-off mi we multiply the size of NLO corrections in our 
NLO fit by mi/nis. The relevant formulae can be easily obtained from standard SU(2) heavy meson 
chiral perturbation theory under the simplification that vector contributions are dropped [05l]. 

We find that SU(3)xSU(3) NLO partially quenched chiral perturbation theory does not 
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Table 4: Fitted parameters from different fits with a valence mass cut am^^g < 0.01. For each fit the LECs 
are quoted at two diffe rent scales A^. (Note: the value of Bq depends on the renormalization scheme like 
the quark masses: to obtain Bq, e.g., in the MS(2GeV) scheme, one has to divide the values quoted here by 
Z^(2GeV). 



(ILs-Ls) Ls (2L6-L4) L4 
SU(3) xSU(3) : aBo = 2.35(16), a/o = 0.0541(40) 

IGeV 5.19(45) -10-4 2.51(99) • lO^^ -4.7(4.2) • lO^^ -6.7(8.0) • lO^^ 

770MeV 2.43(45) • 10-4 8.72(99) • lO^^^ -0.1(4.2) • lO^^ 1.39(80) • lO^^^ 

SU(2) X SU(2) : aBo = 2.414(61), afo = 0.0665(21) 

IGeV 4.64(43) • 10-4 5.16(73) • 10"^ -7.1(6.2) • 10"^ 1.3(1.3) • 10"^ 

770MeV 5.0(4.3) • 10-^ 9.30(73) • lO^^ 3.2(6.2) • 10"^ 3.3(1.3) • lO^^ 



Table 5: Comparison of converted SU(3) x SU(3) fit parameters with those from SU(2) x SU(2) fits. Low 
energy scales [3 4 are defined at A = 139MeV. 







aBo 


afo 


h 


k 


SU(3) 
SU(2) 


X SU(3), conv. 
X SU(2) 


2.457(78) 
2.414(61) 


0.0661(18) 
0.0665(21) 


2.87(28) 
3.13(33) 


4.10(05) 
4.43(14) 



describe our data well up to meson masses comparable to the kaon mass. A good fit can only 
obtained with a cut in bare quark mass of aniavg = {ariix + amy)/2 < 0.01. while the kaon coitc- 
sponds to awjavg — 0.016. The utility of SU(3)xSU(3) ChPT at NLO is questionable; however we 
can quote LECs that fit the data for valence strange and light quark masses obeying the bound 
aniavg < 0.01, and light dynamical quark masses ami < 0.01 but with a fixed dynamical strange 
quark mass of am^ = 0.04. These LECs may well differ from those that would be obtained in the 
unphysical true SU(3) chiral limit. The low energy constants obtained are displayed in table ^ for 
two popular choices of the chiral scale. 

We also directly applied SU(2)xSU(2) to fit the LECs of the effective two flavour theory that 
matches our simulated 2+1 flavour "real" world. These ai^e also displayed in table 0. We fit the NLO 
forms to our data using a mass cut aniavg < 0.01 to obtain a good quality of fit. We perturbatively 
convert our results to the scale independent Ti,, I4 SU(2) LECs in table |5[ Interestingly our 2+ If 
results, both from a perturbatively converted SU(3) fit and from a direct SU(2) fit are broadly 
consistent with each other and with the 2f results of ETMC and CERN gl], This adds 
somewhat to the picture discussed recently by Leutwyler [p3|]. 



5. Quark masses and lattice spacing 

We determine a^\amud and anis from a combination of the S=3 Q. baryon mass and the 
pseudoscalar kaon and pion masses. These quantities then produce the correct Q. mass from 
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Table 6: Determined lattice scale and spacing and unrenormalized quark masses (amx^^^ = am^^'^ +amres)- 

a-^/GeV ^Vfm am^^^^' am^/'' amf^' 

1.729(28) 0.1141(18) -0.001847(58) 0.001300(58) 0.0343(16) 0.0375(16) 



linear extrapolation in the valence mass to anis and to am„^/ in the Ught sea masses. They simul- 
taneously produce consistent kaon and pion masses using our SU(2)xSU(2) chiral extrapolations. 
The scale from the Q. can also be used to predict /;;:, /k, j- and the quark masses Based 
on our preliminary analysis our bare quark masses correspond to renormalised quark masses of 
niud = 3.72(16)MeV and m, = 107.3(4.5)MeV in the MS at 2 GeV. We find = 28.8(4). We 
obtained = 124.1(3.6) MeV, fK = 149.6(3. 6)MeV, which are around 5% lower than their ex- 
perimental values. This is likely an 0{a^) effect, and our ratio j- = 1.205(18). This implies 
\Vus\ = 0.2232(34). Here the eiTors on decay constants ai^e statistical only, and the quark masses 
and Vus contain only partial systematic eiTors. Full systematic eiTors will be estimated in a journal 
paper [p9|]. 

The lattice spacing determined from Mp is somewhat different, being around 1.65 GeV on 
24^ and 1.62GeV on 16^ While on 16^ we relied on vector meson states (which are unstable in 
QCD) and an ad hoc value of 0.495 fm for ro to set the scale, our larger 24^ volume enables fj^ and 
baryon masses to be safely considered. We find broad consistency between decay constants and 
the Omega (which, being composed of three strange quarks, is physically small for a baryon). We 
have found that in the chiral hmit ro/a = 4.13(10), and this suggests we measure a physical value 
of ro ~ 0.47 fm, and disfavour ro = 0.495 fm. Were we to use the pseudoscalar decay constants to 
set the scale ro ~ 0.45 fm. This is a tendency that is consistent with other recent lattice calculations 
[37], and cautions against reliance on vector mesons for precision scale. We are encouraged by 
recent progress that has been made on treating vector meson decay in lattice QCD [44]. 



6. Neutral kaon mixing 

We have updated our paper []T2[] on Bk with 24^ results that have been presented by Cohen and 



Antonio at this conference [^. We use the two wall, operator sink method to gain a spatial volume 
average. We use propagators that are the sum of solutions for periodic and anti-periodic temporal 
boundary conditions to eliminate unwanted round-the-world propagation. This gives exceedingly 
long plateaux on our A'^f = 64 lattice. We have modified our analysis to set the lattice spacing from 
the Q. mass and now use fits assuming only SU(2)xSU(2) chiral symmetry which we consider to 
be more theoretically robust. We have access to lighter masses and more statistically precise data 
and see evidence of curvature in the fixed strange mass chiral extrapolation. Partially quenched 
SU(2)xSU(2) chiral perturbation theory both describes the valence and sea mass dependence well 
in the region of our fit. The smaller volume m/ = 0.02 data point is not included in our fit, but the 
unitary fit curve matches onto this data point reassuringly well. 
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0.58 
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m +m 

X res 

Figure 8: Results for Bp together with the NLO partially quenched SU{2) x SU{2) ChPT fit to the 24^ 
data plotted versus the light valence quark mass m.,. From top to bottom on the left-hand-side, the three 
curves are nii = 0.01, 0.005 and nix respectively. The valence strange quark mass is fixed at its unitary value 
niy = nis = 0.04. While the statistical errors are large, the growing upward curvature in as the sea quark 
mass is increased from 0.005 to 0.01 predicted by ChPT is visible. The values are slightly shifted for 
clarity. 

Non-perturbative renormalisation 2% 



Sea strange mass adjust 1% 

Chiral extrapolation 2% 

Discretisation 4% 

Finite volume 1% 



Table 7: Breakdown of systematic error estimate for our 24^ Bk result. 

The final results follow, with the first error statistical and the second systematic. 

51^(2 GeV) = 0.514(10)(25), (6.1) 
S^(2GeV) = 0.524(10)(28), (6.2) 
Bk = 0.720(13)(37), (6.3) 

The components of the systematic error are shown in table 0. These eiTors are added in quadra- 
ture, and the discretisation systematic is dominant. Simulations in progress with a finer lattice spac- 
ing will directly address this. The two-loop perturbative conversion to MS is cuiTcntly a subleading 
error but will soon become the most important eiTor to address. Finer lattice spacings will only 
yield logarithmic improvement, and a higher order calculation, preferably with non-exceptional 
momenta, is important. A non-perturbative step scaling approach could, of course, even better 
address the convergence of perturbation theory. 
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7. K13 form factor 

James Zanotti presented a status update of our calculation of the semileptonic kaon decay 
form-factor [||], /+(0) which is obtained from the ^ 71 matrix element of the weak vector current 

{7l{p')\V,\K{p))=Mq^)ip^+p'^)+f4q^){p,-p'^) 

This is a promising approach for an accurate determination of Vus, and makes use of standard 



double ratio techniques |45] to measure the deviation of the form factor from unity, giving a very 



small overall error. An example ratio is given below. 



max/ 



Preliminary results using only the niu = 0.03, 0.02, 0.01 data points were presented in [17, 16]. We 
have now added the lightest data point (m„ = 0.005) and finalised our analysis in a full paper |^]. 
The updated analysis includes unified chiral and extrapolations, using a fit form that combines a 
quark mass dependent pole dominance model with the constraints of the AdemoUo-Gatto theorem: 

Jm ,m^,mj^) = 3 . (7.1) 

1 2 

(Mo+Mi(m2+m2))2 

The results from a fit to the large volume (24^) data is presented in Fig. ^ The left and right plots 
show the q^ and quark mass dependencies of Eq. ([0|), respectively. At the physical meson masses, 
we obtain 

/:^'^(0) =0.9644(33) (34) (14) , (7.2) 

which very much favours Leutwyler-Roos [^] results over more recent higher order calculations 
[ p9| , The first error in Eq. ( [7.2| ) is statistical, while the second is an estimate of the systematic 
error due to our choice of ansatz ( [7.1[ ) and the third is the estimate of discretisation errors. The 
PDG quotes |K.v/+(0)| = 0.2169(9) [||]^ so using our result (^, we obtain 



|K.|=0.2247(9)exp(ll)/40) • 

Despite being less mature, K13 form factor appears very competitive with respect to fx/ fn as a 
lattice method for constraining Vus- We anticipate a substantial reduction in error by a recalculation 
with a combination of twisted boundary conditions [|l^] (removing the systematic uncertainty in 
the q^ extrapolation) and stochastic volume averaging for our m„ = 0.005 data point (reducing the 
error in the most important point in the chiral extrapolation). Finally, discretisation effects will be 
addressed before the next lattice conference using our new ensembles with a finer lattice spacing. 

8. Pion and kaon distribution amplitudes 

Chris Sachrajda presented [Q] a calculation of the first and second moments of distribution 
amplitudes of the pion and kaon computed from the following matrix elements 

{K{q)\s{0)y5y{pD^}dmO) = fKiqpiqn{^)K 



more recent analysis finds !V,„/+(0)| = 0.21673(46) 
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Figure 9: We use a unified fit on the 24^ x 64 data to both the and mass dependence of the form factor. 
The data can be usefully displayed in two ways. Left panel contains a pole dominance model interpolation 
of the form factor /o (q^) to q^ =0 having extrapolated to the chiral limit; the data points are adjusted, using 
the fit model, such that if the fit were perfect they would all lie on the fit model. The small remaining scatter 
is an indication of the quality of our unified fit. The right panel shows the chiral extrapolation of fo{q^ = 0) ; 
here the line is our fit model, while the data points are the results of interpolation to = for each ensemble 
and these are consistent. We favour the Leutwyler-Roos prediction and have a smaller error 
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Figure 10: Bare values of {^)k vs the quark mass. The physical region msa—niqa = 0.0375(16) is marked. 



(;r(^)|M(0)757{pD^D^|fif(0)|0) = fniqpiqiiiqv{^'^)K- 



This calculation was performed on our 24^ ensembles and follows on from an earlier work on 16'^ 



[19, hq]. The first moment vanishes in the mass degenerate case and is non-zero for the kaon, 



figure but not for the pion. The second moment has relatively weak mass dependence for both 
kaon and pion, figure 1 1 . These were renormalised using one-loop lattice perturbation theory and 
we obtain the following preliminary results: 



(^)^^(2GeV)= 0.029(2), (,^2^^S(2GeV) = 0.28(3) , ( (2 GeV) = 0.27(2) . (8.1) 
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Figure 11: Bare values of and asa function of the quark mass. 



9. Vector meson decay constants 

These are defined through 

(0 1^2(0)7^^1 (0)|y(p; A)) =fvmve^ 
(0|^-2(0)C7'^V(0)|V(;,;A)) = if^{^i) (£,V-£,V) . 

The vector meson decay consnant fv is well constrained experimentally, but the tensor current 
coupling is useful input that lattice gauge theory can provide to sum rules and other phenomeno- 
logical applications. Chris Sachrajda presented a 2+ If DWF calculation of these transverse decay 
constants, renormalised with Rl-mom NPR using both 16^' and 24^ volumes. 

These results were obtained using only the input strange quark mass of 0.04, rather than the 

fT 

more physically realistic input quai^k mass (i.e. not including wires) of 0.0343. The ratios display 
very weak dependence on m/ and an estimate absorbed this change in sti'ange quark mass to the 
physical point. As the tensor cuiTcnt is scheme and scale dependent the final results were quoted at 
2 GeV in the MS scheme as 

/f (2GeV) Zr(2GeVa) f^^^'-'^a) ^ fp-^-ja) 

f, , 7,, /■bare ' ^ ' rbare ' v • / 

JV Jy Jy 

In the MS scheme with jj. =2 GeV we finally obtain: 

fT fT fT 

^=0.681(20); ^=0.712(11); ^=0.751(9). (9.2) 

Jp JK* J(j> 



10. Nucleon mass and structure 



Takeshi Yamazaki and Shigemi Ohta have presented results [Q] for isovector form factors and 
low moments of structure functions of the nucleon, and related work has also been discussed at 
other conferences [BTI, 52 1. These are performed using our four 24^* ensembles and the correspond- 



ing Edinburgh plot [53| is displayed in figure 12 



Nucleon three point functions have been calculated using a source-sink time separation of 12. 
For our (2.7fm)^ simulation we find, figure 13, that the axial charge gA appears flat except at our 
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Figure 12: Edinburgh plot obtained from various RBC and UKQCD joint ensembles. The DBW2 gauge 
action was used with j3 = 0.764 and 0.72 while the Iwasaki gauge action was used with j3 =2.13 and /3 =2.2. 
The red data points represent the nucleon masses for ot/ G {0. 005, 0.01,0. 02, 0.03}for our 24^ ensembles. In 
the absence of a controlled extrapolation including chiral non-analyticities, our 2+ 1 f results suggest plausible 
agreement with experiment for Nucleon masses, and also suggest reasonable scaling behaviour across several 
couplings and gauge actions. 



lightest datapoint, which is around 15% lower. Similar behavour was seen by RBC, at heavier 
masses on a (1.9fm)^ 2-flavour DWF simulation. We believe this is a finite volume effect with the 
mass threshold determined by the volume. Our m,, = 0.01 data point on our 16^ ensemble does not 



display this effect, but carries very large statistical errors. Figure [13| displays our results against 
trijcL, and is suggestive of mass dependent finite volume effects which appear to scale with mj^L 
and appear for nij^L < 6. We observe similar behaviour if we plot the results with Wilson fermions 
by LHPC/SESAM and QCDSF in this fashion [0, |5|, Hi. Improved statistics for the 16^ 2+lf 
measurements is important to clarify the one (statistically questionable) exception to this picture. 

The lowest 24^ data point is omitted from an extrapolation, and gA = 1.16(6) obtained at the 
physical pion mass. Results were also presented for the vector, axial, induced tensor and induced 
pseudoscalar form factors, some associated couplings (such as g^NN and the induced pseudoscalar 
coupling gp) and corresponding mean squared radii. The momentum fraction, helicity fraction, 
transversity and twist-3 di structure function moments were calculated. 

11. Conclusions 

RBC and UKQCD have exploited the PPARC, Riken, and SciDAC QCDOC machines in Ed- 
inburgh and Brookhaven to simulate dynamical domain wall fermions with realistic sea quark con- 
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Figure 13: Left panel shows our results for g/^ as a function of for several ensembles. Potential finite 
volume effects explain differences between these, and a deviation from experiment. The right panel displays 
the same data as a function of m^^L and the scaling of the deviant points with volume becomes apparent. The 
rather poorly determined blue result for our 16^ ensembles needs more effort to establish whether it confirms 
this picture. 



tent. The analysis of the first lattice spacing a^^ = 1 .73 GeV on a (2.7fm)'' volume is well advanced 
with a broad and rich physics programme presented at this conference. The programme will con- 
tinue to analyse two ensembles on a finer lattice spacing that are currently being generated. The 
physics parameters are very competitive despite the cost of the additional fifth dimension, with 
sea pion masses down to 330 MeV and valence pions down to 240 MeV. Partially quenched chiral 
perturbation theory is exploited in our analysis programme and enabled many more data points to 
be measured within the SU(2)xSU(2) chiral regime. We have obtained results for the low energy 
constants of the chiral effective lagrangian, quark masses, Bk, and Vus from both /k/ fn and from 
K13. We find that very large lattice volumes lUjiL > 6 may be required for non-spectral nucleon 
physics based on suspected finite volume effects in the nucleon axial charge. 
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